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Abstract 

The leading correction from interactions to the transition temperature T c for Bose-Einstein 
condensation can be obtained from a nonperturbative calculation in the critical O(N) scalar field 
theory in 3 dimensions with N = 2. We show that the linear 5 expansion can be applied to this 
problem in such a way that in the large- N limit it converges to the exact analytic result. If the 
principal of minimal sensitivity is used to optimize the convergence rate, the errors seem to decrease 
exponentially with the order in the 5 expansion. For N = 2, we calculate the shift in T c to fourth 
order in 5. The results are consistent with slow convergence to the results of recent lattice Monte 
Carlo calculations. 
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I. INTRODUCTION 



The critical temperature T c for Bose-Einstein condensation in an ideal gas of identical 
bosons was deduced in 1924 by Albert Einstein [Ij]: 

K rp 27T 

fc ^ c - C (| )2 /3 - > W 

where n is the number density of the bosons. An obvious next question is the shift AT C in 
the critical temperature due to the interaction between the bosons. If the potential between 
two bosons is short-ranged, their scattering amplitude at sufficiently low energy is a constant 
proportional to the s-wave scattering length a. The Bose gas with such an interaction has 
been studied extensively since the pioneering work of Lee and Yang in 1957 on bosons 
with a hard-sphere potential. However the shift in T c from such an interaction remained an 
unsolved problem until recently. 

A major breakthrough was made by Baym et al. ||, who showed that the problem could 
be reduced to a nonperturbative calculation in an effective 3-dimensional statistical field 
theory with 0(2) symmetry. This demonstrated conclusively that the leading order shift 
was linear in a: 

ATc -cn^a, (2) 



T, 



c 



where c is a numerical constant. Using a self-consistent calculation of the propagator, they 
estimated the coefficient to be c = 2.9 ||. Baym, Blaizot and Zinn- Justin calculated the 
coefficient analytically in the large- N limit and obtained c = 2.33 [|J. The first correction 
in the 1/N expansion was calculated by Arnold and Tomasik and it reduces c to 1.72 
||. Definitive calculations of the coefficient for the physical case N = 2 were carried out 
by Kashurnikov, Prokof'ev, and Svistunov || and by Arnold and Moore [|7| by applying 
Monte Carlo methods to a lattice formulation of the effective field theory with the result 
c = 1.32 ± 0.02. The second order correction to AT C /T C proportional to (an 1 / 3 ) 2 has since 
been calculated by Arnold, Moore and Tomasik ||. Before the recent breakthrough, there 



had been many previous attempts to calculate the leading shift in T c [^j, [T(], 11, 12 1. The 
results for the coefficient c in (^) ranged from —1 to +4.7. The definitive solution of this 
problem provides a way to assess the reliability of the various calculational methods. 

One particularly simple method that can be applied to this problem is the linear 5 
expansion ||13|| , also known as optimized perturbation theory |14| or variational perturbation 
theory [THj. With this method, an arbitrary parameter m is introduced into the theory and 
calculations are carried out using perturbation theory in a formal expansion parameter 5 
whose value is 1. The rate of convergence can be improved by adjusting the parameter m 
at each order in 5 by some criterion such as the principal of minimal sensitivity. It was 



first applied to the problem of calculating AT C by de Souza Cruz, Pinto and Ramos |TB| . At 
second, third, and fourth orders in 5, the predictions are c = 3.06, 2.45, and 1.48, respectively 



17|| , which seem to be converging to the lattice Monte Carlo result. 

In this paper, we revisit the application of the linear 5 expansion to the calculation of 
the shift in T c . In Section II, we review the reduction of the calculation of AT C to a problem 
in a 3-dimensional statistical field theory. In section III, we introduce the linear 5 expansion 



and review the results of Refs. [|T^]. We present two alternative precriptions for applying 



this method to the shift in T c . In Section IV, we show that in the large-iV limit, these 
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alternative prescriptions converge to the analytic result of Ref. || for some range of m. We 
study the improvement in the convergence rate that is obtained by using the principal of 
minimal sensitivity to choose a value of m that depends on the order in 5. We show that 
there is a family of solutions to this equation that seems to improve the convergence rate 
from a power law in the order of the S expansion to exponential. In Section V, we calculate 
the result for iV = 2 at 3 rd and 4 th orders in 5. In Section VI, we compare our results 
with the previous calculations using the linear 5 expansion of Refs. [ IB , 17 ] . Our results are 
consistent with slow convergence to the Monte Carlo results of Refs. 0, [7J. 



II. EFFECTIVE THEORY 

We consider a system of identical bosons of mass m at a temperature T and number 
density n close to the critical temperature T c for Bose-Einstein condensation. If the effects of 
interactions are small, Einstein's formula (H) implies that the interparticle spacing n _1//3 and 
the thermal wavelength Ay = (2nh 2 /mkBT) 1 / 2 are comparable. We assume that they are 
both large compared to the range R of the two-body potential. In this case, the interaction 
between the bosons can be characterized entirely by the s-wave scattering length a, which we 
assume to have magnitude of order R. Thus we require a <C n^ 1 ^ 3 , At- This system can be 
described by a second-quantized Schrodinger equation. The corresponding imaginary-time 
Lagrangian is 

c = - ^*v 2 v - i*w + (3) 

or 2m m 
where /x is the chemical potential associated with the particle number density. In this paper, 
we generally use units such that H = kg = 1. Treating the system at finite temperature 
using the imaginary time formalism, the field ?/>(x, t) can be decomposed into Fourier modes 
ip n (x) with Matsubara frequencies ui n = 2nnT. The zero frequency mode ^o( x ) is defined 
by 

V>o(x) = T j 1,T dr^t). (4) 



o 

If the temperature is near the critical temperature T c , the chemical potential \x is small in 
magnitude compared to T. For momenta small compared to (mT) 1//2 , the nonzero frequency 
modes have Euclidean energies of order T, while the zero frequency modes have Euclidean 
energies of order Thus at large distances much greater than Ay, the nonzero modes 
decouple from the dynamics, leaving an effective theory of only the zero modes [Q]. A naive 
derivation of the effective lagrangian is obtained simply by substituting ^(x.,t) — > ^o( x ) i n 
(0). The action then reduces to 

fJ T dr J d 3 x £[V>(x,t)] — > J d 3 x £ eff [V>o(x)], (5) 
where the effective lagrangian is 

£ eff = -^^oVVo - £Wb + ^(«o) 2 . (6) 

It is convenient to replace the complex field i/jq by the 2-component real field <fi = (^1,^2) 
defined by 

^ o (x) = (mT) 1 / 2 (0 1 (x) + z0 2 (x)). (7) 
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The effective lagrangian density for the dimensionally reduced theory can then be written 

as 

£ cff = -\$. V 2 + \4 2 + (0 2 ) 2 , (8) 

where the parameters are 

r = —2m\i, (9) 
u = 48vr amT. (10) 

This 0(2) field theory can be generalized to an O(N) field theory by replacing by an 
iV-component real field <fi = (0i, 0at). 

This naive derivation of the effective lagrangian is supported by a more sophisticated 
analysis within the effective field theory framework M. This analysis reveals that there are 
also higher dimension operators in the effective lagrangian (H), such as ((f) 2 ) 3 and V 2 0- V 2 0, 
but they are not needed to calculate the shift in T c to leading order in an 1 / 3 . It also reveals 
that there are corrections to the coefficients r and u that can be calculated perturbatively 
as expansions in o/Ar- Again, they are not needed to calculate AT C to leading order in 
an 1 / 3 . Finally, since the scale T has been removed from the problem, the effective field 
theory defined by the effective lagrangian (§) requires an ultraviolet cutoff. Fortunately, 
the leading order shift in T c is dominated by the infrared region of momentum space and is 
insensitive to the ultraviolet cutoff. 

From Einstein's formula ([I]), the leading order shift AT C in the critical temperature at 
fixed number density is related to the leading order shift An c in the critical density at fixed 
temperature by the simple relation 

AT C 2An c . . 
- = -. 11 

T c 3 n c 1 ; 

The number density in the original theory is n = , where the angular brackets indicate 

the expectation value in the state corresponding to a homogeneous gas. The number density 
in the effective theory, to the accuracy required to calculate the shift in T c to leading order 
in an 1 / 3 , is 

n = n + mT (<f) 2 } , (12) 

where n is the short-distance contribution to the number density, which depends on the 
ultraviolet cutoff. The angular brackets indicate the statistical average of the operator <fi 2 . 
It can be expressed as 

(0 2 ) = 2/ L 2 + r + S(p)] _1 , (13) 

where is the self-energy of the <fr field. The integral over p is divergent and requires an 
ultraviolet cutoff. The number density for an ideal gas is given by the same expression with 
a = 0, which implies S(p) = 0: 

n\ a=0 = n + 2mT f p 2 + r . (14) 

Again, the integral over p is divergent and requires an ultraviolet cutoff. To the accuracy 
required to calculate the shift in T c to leading order in an 1 / 3 , the constant n is the same 
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as in (|12"D . The critical point is determined by the condition that the correlation length 
associated with the propagator in (O) is infinite: 

r + S(0)=0. (15) 

This condition reduces to r = if a = 0. The change An c in the critical density at fixed 
temperature T due to interactions is the difference between (|T^) with r = — S(0) and ( |H]) 
with r = 0. The expression flTT| ) for the fractional shift in T c can therefore be written as 

AT C 2 ttiT 

-zr- = A > 16 

T r 3 n ' v ; 



where A is a quantity that can be calculated within the effective field theory: 



N 



p 2 + E(p) - £(0) 



1 V^ J 



(17) 



The integral over p is convergent and therefore no longer requires an ultraviolet cutoff. We 
have replaced the prefactor 2 in ([T7|) by iV to generalize it to the O(N) case. Comparing 
( PH) with (0) and using ( JlQ| ) and ([!]), we find that the coefficient c is 

3 U2j (-n/487T 2 )- 1 j 

The effective theory defined by the lagrangian (|3|) is superrenormalizable. Although the 
theory has ultraviolet divergences, A is dominated by infrared physics and it is insensitive to 
the ultraviolet cutoff. The relevant length scales are therefore those set by the parameters r 
and u of the effective field theory. The critical point can be reached by tuning the parameter 
r. At the critical point, there is only one relevant length scale and it is set by the parameter 
u. Since A has dimensions of length, it must be proportional to u by simple dimensional 
analysis. The determination of the coefficient of u requires a nonperturbative calculation. 
Since u is linear in a as a — > 0, this demonstrates that the leading term in the shift in T c is 
linear in a. 

The calculation of the first order shift AT C in the critical temperature has been reduced 
to the calculation of the quantity A in (|17D for N = 2. The first calculation of this quantity 
within a systematically improvable framework was a calculation to leading order in the 1/N 
expansion by Baym, Blaizot, and Zinn- Justin jp: 

Nu 

A = --— (large N) . (19) 

The next-to-leading order (NLO) correction in the 1/N expansion was calculated by Arnold 
and Tomasik || . The first two terms in the 1/N expansion are 

/ 0.527 \ / Nu \ 
A =( 1 - — + -)(-96^)' < 2 °> 

For the physical case iV = 2, A has been calculated nonperturbatively by Kashurnikov, 
Prokof'ev, and Svistunov and by Arnold and Moore [0] using Monte Carlo simulations 
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of a lattice formulation of the 0(2) field theory. The two results are consistent, but Ref. |7j 
quotes a smaller error: 



A = (0.568 ± 0.008) [-^) (N = 2) . (21) 

If we set N = 2 in the 1/N expansion in fl20|) , the coefficient of — u/487r 2 is 1 at leading 
order and 0.74 at next-to-leading order. These coefficients seem to be approaching the lattice 
Monte Carlo result 0.57. 



III. LINEAR 6 EXPANSION 



Any attempt to calculate AT C using ordinary perturbation theory in u is doomed to 
failure. In the symmetric phase where r > 0, perturbation theory is an expansion in powers 
of u/r 1 / 2 . If we naively let r — > to find the critical behavior, the perturbative corrections 
are infrared divergent, with the degree of divergence becoming more severe at each successive 
order. It is useful to think of perturbation theory as an expansion in powers of u£, where 
£ is the correlation length. Perturbation theory breaks down in the critical region, because 
w£ is of order 1 and every order of perturbation theory gives a leading order contribution. 

The linear 5 expansion (LDE) is a particularly simple method for obtaining nonpertur- 
bative results using perturbative techniques. Calculations are carried out as in ordinary 
perturbation theory by evaluating a finite number of Feynman diagrams. The name "linear 
5 expansion" comes from the fact that it can be defined by a lagrangian whose coefficients 
are linear in a formal expansion parameter 5. If £ is the lagrangian for the system of interest, 
the lagrangian that generates the LDE has the form 

£ s = (l-5)£ + 5£, (22) 

where £q is the lagrangian for an exactly solvable theory, The lagrangian £$ interpolates 
between £ Q when 5 = and £ when 5 = 1. The linear 5 expansion is much easier to apply 
than the nonlinear 5 expansion in which the lagrangian that interpolates between £ and £ 
involves an operator that is a nonlinear function of 5 [^3[ . 



The lagrangian £$ for the exactly solvable theory typically involves an arbitrary param- 
eter m. At 5 = 1, the lagrangian fl22|) is independent of m. If a perturbation series in 5 
converges at 5 = 1, it must converge to a value that is independent of m. However at any 
finite order in the 5 expansion, results depend on m. As m varies over its physical range, 
the prediction for the observable often extends out to ±oo. Thus some prescription for m 
is required to obtain a definite prediction. A simple prescription for m is the principle of 
minimal sensitivity (PMS) that the derivative with respect to m should vanish. The LDE 
supplemented by the PMS prescription is often referred to as optimized perturbation theory 
|T4| or variational perturbation theory jnj|. A prescription such as PMS may improve the 



convergence rate of the LDE, because it allows the variable m to change with the order n 
of the LDE. 

The convergence properties of the LDE have been studied extensively for the quantum 
mechanics problem of the anharmonic oscillator. The hamiltonian is 

H = \p 2 + \uj 2 q 2 + gq A - (23) 
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The LDE is generated by the hamiltonian H$ = (1 — 8)Hq + 8 H, which interpolates between 
Hq at 8 = and H at 8 = 1. It can be expressed as H = Hq + -ffmt, where 

#o = ~p 2 + \m\\ (24) 

#int = ~^(m 2 - u 2 )q 2 + 8gq\ (25) 

The LDE has been proven to converge for appropriate order-dependent choices of the variable 
m that include the PMS criterion as a special case. Duncan and Jones showed that the 
finite temperature partition function converges exponentially, with the errors at n th order 
decreasing as exp(— bTn 2 ^ 3 / g 1 ^ 3 ) where b is a numerical constant |]I9 |. Guida, Konishi, and 



Suzuki proved that the energy eigenvalues E n converge uniformly in g as n — > oo [[2(J . Janke 
and Kleinert showed that the leading coefficient in the strong-coupling expansion for the 
ground state energy E converges exponentially, with the errors decreasing like exp(— cn 1 ^ 3 ) 
where c is a numerical constant |21| . 



Belief, Garcia, and Neveu have made numerical studies of the LDE optimized by the 



PMS criterion for the anharmonic oscillator |[22|| . The PMS criterion is a polynomial equation 



in the variable m, so its solutions are in general complex-valued. The authors of Ref. [22 
found that the solutions at different orders n in the LDE can be organized into distinct 
families. The most accurate predictions for the ground state energy E and for the quantity 
{q 2 ) are obtained by using families of complex solutions. The resulting values of E and {q 2 ) 
are of course complex-valued, but the imaginary parts are small. The authors of Ref. f2"2fl 
suggested that the real parts could be taken as the predictions of the LDE and the imaginary 
parts could be used as an indication of the error. 

The anharmonic oscillator can be regarded as a Euclidean field theory with a single 
real-valued field in 1 space dimension. The statistical field theory defined by the lagrangian 
(§) is a generalization to a multicomponent field in 3 space dimensions. The most serious 
obstacles to generalizing the convergence proofs for the anharmonic oscillator to this more 
complicated problem come from the infrared and ultraviolet regions of momentum space. 
It is reasonable to expect the convergence behavior to be similar if appropriate infrared 
and ultraviolet cutoffs are imposed on the field theory. The quantity A defined in ([H]) is 
insensitive to the ultraviolet region, so we do not expect any complications from taking the 
ultraviolet momentum cutoff to oo. However, A is very sensitive to the infrared region, so 
convergence in the presence of an infrared cutoff gives no information about the behavior of 
the LDE in the limit when the infrared momentum cutoff goes to 0. 

We proceed to apply the LDE to the O(N) statistical field theory defined by the la- 
grangian C c s in (|8]). Our choice for the solvable field theory is the free field theory with 
mass m. The lagrangian defined by setting C = £ fr in ( |2"2"D can be written £ = £ + £mt, 
where 

£ = _i^.v 2 ^+im 2 ^ 2 . (26) 

C mt = S -(r-m*)$*+ S ^(P) 2 . (27) 

Calculations are carried out by using 5 as a formal expansion parameter, expanding to a 
given order in 5, and then setting 5 — 1. Thus the Feynman diagrams generated by the LDE 
will involve a propagator with mass parameter m 2 , a 4-point vertex proportional to 5 u, and 
a mass counterterm proportional to 8(r — m 2 ). 
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In order to apply the LDE to the shift in T c , we need a prescription for generalizing 
the quantity A defined in (|TTD to the field theory whose lagrangian is £5 = Co + C\ nt . The 
prescription must have a well-defined expansion in powers of 5, and it must reduce to ( |I~7D 
when 5=1. One possibility, which we will refer to below as prescription II, is to use the 
expression (|TT|), where is the self-energy for the field theory defined by £5. Any such 
prescription defines A(u, m, 5) as a function of three variables. The n th order approximation 
in the LDE is then obtained by truncating the expansion in powers of 5 at n th order to 
obtain a function A^(u, m, 5) and then setting 5=1: 

A « A {n) {u,m,5 = 1). (28) 

This sequence of approximations may not converge as n — > 00, but if it does converge, we 
expect it to converge to the correct value A. 

At any finite order n, the approximation (pgj) depends on m. We will see that as m 
varies over its physical range from to +00, the range of the predictions for A extends out 
to +00 for n even and out to —00 for n odd. Thus some prescription for m is required to 
obtain a definite prediction. One prescription is the principle of minimal sensitivity (PMS): 

PMS : ^-A^ n \u,m,5 = 1) = 0. (29) 
dm 

After setting 5 = 1, A*™-* is a function of u and m only. By dimensional analysis, a solution 
m to the PMS criterion must be proportional to u and the value of A*™) at such a point must 
also be proportional to u. An alternative prescription is the fastest apparent convergence 
(FAC) criterion: 

FAC: A( n »(ti ) m,i = l)-A<"- 1 '(i 1) m^=l) = 0. (30) 

By allowing the variable m to change with the order n, prescriptions such as PMS or FAC 
may improve the convergence rate of the LDE. 

The LDE was first applied to the problem of calculating the shift in T c by de Souza 



Cruz, Pinto and Ramos [16]. Their prescription for the quantity A defined in (|T7|) is 



Al = N J p [p 2 + m2 ( 1 " 5 ) + S G°) '-(V) 1 > ( 31 ) 

where S(p) is the self-energy for the field theory whose lagrangian is £,&. This prescription 
agrees with (p~7|) at 5 = 1. Note that in addition to the explicit dependence on 5 in (0), there 
is also implicit dependence on 5 through S(p), which can be calculated as a power series in 5. 
The subscript I on A in ([31]) distinguishes this prescription from alternative prescriptions for 
A that we will label II and III. In the Feynman diagrams for evaluating S(p) in the LDE, all 
the propagators have mass m. Thus when the prescription I for A is expanded in powers of 
5, all the integrals involve only propagators of mass m, with the exception of the subtraction 
term in (|3~i"D , which involves a massless propagator. The quantity A was calculated to 2 nd 
order in 5 in Ref. [[RJ and to 3 rd and 4 th order in Ref. |L7|]. After truncating the expansion 
in powers of 5 and setting 5 = 1, they determined m by applying the PMS prescription 
(p5p. At n th order in the 5 expansion, this reduces to a polynomial equation in m with n 
roots, most of which are complex- valued. For a particular choice of the roots at 2 nd , 3 rd , 
and 4 th order in 5, the real parts of A/(— u/48tt 2 ) are 1.31, 1.05, and 0.68, respectively. 
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This sequence seems to be approaching the lattice Monte Carlo result 0.57. In Ref. ||17|| , the 
authors also applied the LDE to the second order shift in T c with encouraging results. 

Since the result for A is known analytically in the large N limit, one would also like to 
compare the predictions of the LDE to this analytic result. Unfortunately, as we shall see, 



the prescription I of Ref. |T7fl is not well-behaved in this limit. We therefore present two 



alternative prescriptions for A that are well-behaved in the large- N limit 

(p 

(p 2 + m 2 (1 - 5) + E(p) - S(0)) _i - (f + m 2 (1 - 6) 



An = iv|[(p 2 + S(p)-S(0)) _1 <- 2 



A m = N 



(32) 
(33) 



These both agree with (|17]) at 5 = 1. With prescription III, the LDE for A generates 
integrals involving only propagators of mass m. With prescription II, the integrals involve 
both propagators of mass m and massless propagators. 

Our prescription III in ( |33"D is closely related to the prescription I in Ref. [|I7| . The 
difference between them is given by a simple integral that is evaluated in ( |A.22 ): 



N 



Ai - Am = -—mVT^5. (34) 

47T 



Because of the additional term in Aj proportional to Ny/1 — 5, the limit 5 — > 1 does not 
commute with the large- N limit. At n th order in the 5 expansion, the terms that dominate 
as iV — > oo are the Nm term and the highest order term in the 5 expansion which is 
proportional to N n u n /m n ~ 1 . The PMS criterion gives a value of m that scales like N 1 ~ 1 ^ n u 
and the corresponding value of A scales like A r2_1//n M. The correct result is that A scales like 
Nu in the large N limit. Thus prescription I has the wrong behavior in the large- N limit 
at any order in the LDE. 



IV. LARGE N LIMIT 

In this section, we study the behavior of the linear S expansion (LDE) in the large- A" 
limit defined by N — > oo, u — > with Nu fixed. We calculate A to all orders in 8 using 
prescriptions II and III. We show that the LDE converges to the correct analytic result 
for m/Nu greater than some minimum value, with errors that seem to decrease roughly 
like n -1 / 3 . We also study the improvement in the convergence obtained by using the PMS 
criterion for m. If we use the family of real solutions of the PMS criterion, the errors in A 
seem to decrease like n" 1//2 . If we use a particular family of complex solutions, the errors 
seem to decrease exponentially with n. 

The quantity A defined in (|I7]) has the diagrammatic representation shown in Fig. [I], 
where the line with the shaded blob represents the propagator [p 2 + S(p) — S(O)] -1 . In the 
large-A limit, the leading contribution to the subtracted self-energy S(p) — S(0) comes from 
the series of bubble-chain diagrams whose 1 st and 4 th members are shown in Figs. and || 
Since these diagrams are of order 1/N, the leading contribution at large N is obtained by 
expanding the expression (|i~7D for A to first order in S(p) — S(0). For the prescription II 
given in (|32"D, the resulting expression is 

Au—NfW-m. (35 ) 
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FIG. 1: Diagram for (cf) 2 ) at the critical point. The black blob represents the cf) 2 vertex, while the 
shaded blob represents the complete propagator with subtracted self-energy S(p) — 2(0). 




FIG. 2: The diagram that contributes to — S(0) at order 8 2 



The corresponding expression for the prescription III is obtained by replacing p 2 — > p 2 + 
m 2 (l — 5) in the denominator of (|35|). Although is ultraviolet divergent, the difference 
— S(0) is finite. At any finite order in the LDE, — S(0) scales as p 2 at small p 
and as Inp for large p, so the integral in fl35|) is convergent. In addition to the diagrams 
for £(]?) generated by the interaction term 5m0 4 in (|27l), we must also take into account 
insertions of 5r = —5 S(0) and insertions of 5 m 2 . The effect of the 5r insertions is to replace 
each 1PI self-energy subdiagram in by the subdiagram with the zero-momentum limit 
subtracted. The effect of the Sm 2 insertions is to replace m 2 by m 2 (l — 5). 

We first consider the prescription II for A. The expression for the large- N diagram for 
A with n + 1 loops is given in ( |A.23| ). The integral I n is a constant multiplied by m~^ n ~ l \ 
After taking into account the 5 m 2 insertions, the sum of the large- N diagrams is 



A n = 2 £ (-8-u 

n=2 v u 



In (1 - S) 



-(n-l)/2 



(36) 



The integrals I n are expressed as 1-dimensional integrals in ( p.23|) . Inserting these expres- 




FIG. 3: The 4 in the series of diagrams for S(p) that survive in the large-^ limit. 
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sions into (|3~6"D, we obtain 



1 \ n—l 

8 l \ f°° y 



2 



A n = 8 -mT j- / dy »- , [A (y)] , (37) 

where /x is the dimensionless mass parameter 

" = -Tsr < 38) 

and A(y) is the function 

2 y 

A(y) = - arctan — . (39) 

The prediction for A at n th order in the 8 expansion is obtained by expanding the expression 
(|37|) as a power series in 8, truncating after order 8 n , and then setting 6=1. 

It is easy to show that if the LDE for (|3"7|) converges, it converges to the correct analytic 



result (0). After interchanging the order of the sum and the integral, the sum can be 
evaluated and we obtain 

r N f°° , y 2 § A (y) / a n\ 

An = S—^u / dy » - — ^ — . 40 

24tt 3 Jo y (yi + l) 2 VT^5 + 5 A(y) V ; 

Upon taking the limit 5 — *■ 1, all dependence on the mass parameter \x disappears: 

A n = -5 -u / dy ? . (41) 

24vr 3 io (y 2 + l) 2 



After evaluating the integral over y, our final result is (|T^) . Note that no prescription for 
/i was required to get this result. Thus, if the 5 expansion converges at some value of /i, it 
should converge to fllPP- 

The manipulations that were used to show the convergence to (O) involved several 
interchanges of limits. First, we summed up some terms in the 5 expansion to all orders 
to obtain the sum over diagrams in (^). Next, we interchanged the order of the sum over 
diagrams and the integral over the last momentum for each diagram. Next we took the limit 
8 —*■ 1. Finally we evaluated the last momentum integral. It is difficult to translate the 
conditions for the validity of each of these steps into a statement about the convergence of 
the LDE. However, the convergence can be easily studied numerically. In Fig. [|, we show 
A as a function of the dimensionless variational parameter /i for several orders in the 8 
expansion: n = 3 and n = 2 3 ' + 3, j = 0, ...,5. The horizontal line is the analytic result 
(19). The results are consistent with convergence to flT9|) for all /x greater than a critical 



value n c which is close to 0.7. Our best estimate for the endpoint of the convergence region 
is jji c = 0.706 ± 0.003. If /x < /x c , A seems to diverge to — oo for n even and +oo for n odd. 
For any fixed value of /x > /x c , the convergence with n is very slow. In column 1 of Table |, 
we give the values of A at /x = 1.039, which is the location of the minimum of the n = 3 
curve in Fig. [|. We define the fractional error e n by 



| Re A - (-Nu/96 



7T 



2\ 



-AV96tt 2 ) ' ^ 
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FIG. 4: A/(— Nu/96tt 2 ) in the large- N limit as a function of [i at n th order in the LDE for n = 3, 
4, 5, 7, 11, 19, 35. The curves for n = 7, 11, and 19 appear in order between those labelled n = 5 
and 35. 

(Taking the real part of A is superfluous for real values of /i, but it will be important later 
when we consider complex values of /j,.) In Fig. [5], we show a log-log plot of e n as a function 
of n. The points lie close to a straight line, indicating that the errors decrease like a power 
of n. The dotted line is e n = 0.70n -0 ' 37 , which goes through the points for n = 19 and 35. 
Thus the errors seems to decrease roughly like n -1 / 3 . 

The rate of convergence can be improved by changing the value of fi at each order in 
the LDE. One possibility is to use the PMS criterion (|2~9"D to determine fi. At n th order in 
the 5 expansion, this condition reduces to a polynomial in fi of order n — 2. For n even, 
there are no real roots. For n odd, there is always one real root. These roots correspond to 
the minima of the curves in Fig. |j. The values of A/(—Nu/9Qtt 2 ) are given in column 2 of 
Table I. There is slight improvement compared to column 1 in the convergence toward the 
correct value 1. The fractional errors defined in (f42"[) are shown as a function of n in Fig. |5|. 
The points lie close to a straight line, indicating that the errors decrease like a power of n. 
The dashed line is e n = 0.78n -0 ' 46 , which goes through the points for n = 19 and 35. Thus 
the errors seem to decrease roughly as rT x l 2 . 

The PMS criterion has either zero or one real solution, depending on whether n is even or 
odd. However there are always n — 2 complex- valued solutions. In the case of the anharmonic 
oscillator, there are families of complex solutions with much better convergence properties 



than families consisting of purely real solutions ||22j| . In our problem, at any odd order 
n, the real solution always has the largest value of Re A, which is the value farthest from 
the correct result (^). Thus this family of solutions gives the slowest possible convergence 
rate. However there is a strong anticorrelation between the errors in ReA and ImA. This is 
illustrated in Fig. ^|, which is a scatter plot of |ImA| vs. ReA for the solutions to the PMS 
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FIG. 5: Log-log plot of the fractional error e n as a function of the order n in the LDE. The squares, 
diamonds, and circles are for fi = 1.039, the real solution to the PMS criterion, and the imaginary 
solution with maximal |ImA|, respectively. The lines are simple curves that pass through the last 
2 points. 

criterion for n = 35. The solutions that give the most accurate values for Re A are those 
with the largest values for |ImA|. Thus we can define a nearly optimal family of solutions 
by choosing those with the maximal values of |ImA|. The values of A for this family of 
solutions are given in column 3 of Table I. Note that ImA for these solutions shows no sign 
of converging to 0, and in fact iseems to be converging to a value close to 0.4 ReA. However 
the improvement in the convergence of ReA compared to the real values of A in column 2 
is evident. A log- log plot of the fractional errors in ReA is shown in Fig. The downward 
curvature of the points indicates that the errors decrease faster than any power of n. The 
solid line is e n = 0.32(0.93)™, which goes through the points for n = 19 and 35. Up to the 
order to which we have calculated, the errors are decreasing faster than this exponential. 

We obtain similar results if we use the prescription III for A given in (|3"3"D. The sum of 
the large- N diagrams is the same as in (|36|), except that the integrals I n are replaced by J n . 



The integrals J n are expressed as 1-dimensional integrals in (p.24| ). Interchanging the order 



of the sum and the integral and then evaluating the sum, we obtain 
a x N /■*> y 2 5 A{y) 

Am = — o -u / ay— — — r-r . . 43 

16tt 3 Jo y (y^ + l)(y^ + A)^/T^Sfi + 5A(y) 1 ; 

Taking the limit 5 — > 1, all dependence on \x disappears and after evaluating the integral 
we again obtain (|i~9"|). Thus, if the 5 expansion converges, it should converge to the correct 
result. 

The convergence can be again studied numerically. We find that it seems to converge 



to (plj) for all n > n c and it seems to diverge for fi < fj, c , with fj, c w 0.7. In column 4 of 
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Re A /(-Nu/96% ) 

FIG. 6: Scatter plot of |ImA| vs. ReA for the solutions fi of the PMS criterion at 35 th order in the 
LDE. 



Table |, we give the values of A for various orders n at fi = 0.947, which is the location 
of the minimum of A for n — 3. The convergence at a fixed value of /i is very slow. As 
in the case of prescription II, the errors seem to decrease roughly as n~ 1//3 . More rapid 
convergence can again be obtained by using the PMS criterion for /i. The PMS criterion 
has no real solutions for n even and 1 real solution for n odd. The values of A for the real 
solution at various orders n are given in column 5 of Table I. This is the family of solutions 
with the slowest convergence rate for ReA. The errors seem to decrease roughly as n -1 / 2 . 
A family of solutions with a nearly optimal convergence rate for ReA consists of those with 
the maximal value of | Im A | . The values of A for this family are shown in column 6 of Table 
I. The errors in ReA are numerically larger than those with prescription II, but they also 
seem to decrease exponentially with n. 



V. FINITE N 



In this section, we use the linear S expansion (LDE) to calculate A for finite N. We 
calculate A to 4 th order in 5 using both of the prescriptions II and III. Setting N = 2 and 
using the PMS criterion, we obtain the predictions relevant to the shifts in T c for Bose- 
Einstein condensation. 

The first term in the LDE for A is second order in 5. It is obtained by calculating the 
3-loop diagram involving the 2-loop self-energy diagram in Fig. 0. The contribution to A 
of order 5 3 is obtained by calculating the 4-loop diagram involving the 3-loop propagator 
correction in Fig. [7|. The contribution of order 5 4 is obtained by calculating the 5-loop 
diagram involving the 4-loop propagator corrections in Fig. The diagrams are reduced to 
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a 



b 



c 



FIG. 7: The diagrams that contribute to S(p) — S(0) at order 8 3 . 

momentum integrals in Appendix [A] and the momentum integrals are evaluated in Appendix 
O. The complete expression for A for the O(N) theory through order 5 4 is 



A n = 5 [ 



,N(N + 2) 



u 2 I 2 1 + -5+-5* 



3 N(N + 2)(N + v, :i 
108 



648 



(N 2 + 6N + 20)/, 



+2(5N + 22) {h b + J 4c ) + 2{N + 2) (3/ 4d + he) 



(44) 



Inserting the values for the integrals from Appendix B, truncating at successively higher 
orders in 5 and then setting 5 = 1, we obtain the results at 2 nd , 3 rd and 4 th orders in the 
LDE: 



A(2) 
^11 



A 



(3) 



.(4) 
^11 



Nu 
~96vr 2 

Nu 
~96tt 2 

Nu 
'96tt 2 



F ■ 

'1 N + 8 1 > 

'51 - 1.03886^ 1 
4 /i JV + 2 /i 2 



+ 



0.430513 



TV 2 + 6N + 20 
(iV + 2) 2 



+ 1.09 



5N + 22 
(iV + 2) 2 



- 0.187713 



iV + 2 



where /i is the dimensionless mass parameter 



487rm 
(iV + 2)w' 



(45) 
(46) 



(48) 



In the large- N limit, this reduces to the mass parameter introduced in (]38|) . The accuracy 
of the 4 th order result in (£|7|) is limited by the accuracy to which we have evaluated the 
integrals I#> and I± c . In Fig. || we show A for N = 2 as a function of \i for n = 2, 3, 4. The 
horizontal lines are the edges of the upper and lower error bars of the lattice Monte Carlo 
result given in (|21|). 
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FIG. 8: Four-loop diagrams that contribute to — S(0) at order 5 4 



At order 5 2 , the PMS criterion has no solutions, because fl45[) is a monotonic function of 
fi. At order <5 3 , the PMS criterion has a single real solution /i (3) = 1.03886(iV + 8)/(N + 2) 
and the resulting value for A is 



A§ } = 0.481296 



N + 2 

ivTs 



Nu 
96V 



(49) 



At 4 th order in 5, A is a monotonic function for real-valued a, so the PMS criterion has 
no real solutions for a. It does however have a conjugate pair of complex solutions whose 
values at N = 2 are u = 2.0772 ± 1.773z. The corresponding values of A are complex, but 
we can take the real part as the prediction for A at 4 th order in the LDE: 



ReA 



(4) 
11 



0.214 



(50) 



For N = 2, the 3 rd and 4 th order approximations in (f49|) and (|50| ) differ from the lattice 
Monte Carlo result of ([H]) by about 66% and 62%, respectively. These predictions are not 
as accurate as in the large N limit, where the errors in the 3 rd and 4 th order approximations 
are about 52% and 44%, respectively. Furthermore, the percentage decrease in the error in 
going from 3 rd order to 4 th order is only about half as large for iV = 2 as it is for large N. 

With the prescription III, the expression for A to 4 th order in 5 is the same as (|44|). 
except that the integrals I n are replaced by J n . Inserting the value for the integrals J n from 
Appendix B, the results at 2 nd , 3 rd and 4 th orders in the LDE are 



.(2) 
^11 



A (3) 

A 



in 

(4) 



Nu 
~96tt 2 

Nu 
'96tt 2 

Nu 



in 



96vr 2 
+ 



0.575364- 

0.863046- - 
v I 1 

1.07881 - - 0.817603 



0.408802 



N + 8 1 
iVT2^ 2 " < 
A^ + 8 1 
'N + 2~^ 



(51) 
(52) 



N 2 + 6N + 20 5N + 22 1 

0.316944 — — ^ — + 0.81 ,,. _ - 0.149896 



(N + 2Y 



(N + 2) 



N + 2 



(53) 



The accuracy of the 4 th order result in (53) is limited by the accuracy to which we have 
evaluated the integrals and J^ c . At order 5 2 , the PMS criterion has no solution. At 
order 5 3 , it has one real solution for fi. At 4 th order in 5, it has a conjugate pair of complex 
solutions. The corresponding values of A are given in Table 0. The values of A at 3 rd order 
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FIG. 9: A/(— u/A8tt 2 ) for N = 2 as a function of /i at n th order in the LDE using prescriptions I 
(upper 3 curves on the right) and II (lower 3 curves on the right). The dotted, dashed, and solid 
curves are for n = 2, 3, and 4, respectively. 

and of Re A at 4 th order are close to the corresponding values from prescription II, but a 
little farther from the lattice Monte Carlo result. 

The results of the LDE for prescription I can be obtained from those for prescription III 
by using (|34]): 

In Fig. ^ we show A*™) as a function of /i for n = 2, 3, and 4. The values of A^ n ^ at the 
solutions to the PMS criterion are given in column 1 of Table |J. At orders 3 and 4, the 
PMS criterion has a conjugate pair of complex solutions. As the predictions of the LDE, 
the authors of Ref. |17| chose the positive value for A^ 2 ^ and the real parts of the complex 
values for and AT 4 ). These predictions differ from the lattice Monte Carlo results by 
about 131%, 85%, and 20%, respectively. 



VI. CONCLUSIONS 

The calculation of the shift in T c to leading order in an 1//3 can be reduced to a non- 
perturbative calculation of a quantity A in a 3-dimensional 0(2) field theory at its critical 
point. This quantity can be generalized to an O(N) field theory and it can be calculated 
systematically using the 1/N expansion. We have shown that the linear 5 expansion (LDE) 
can be applied to the calculation of A in such a way that the large N and 8 — > 1 limits 
commute. The prescriptions II and III defined in (|32|) and (|33"D have this property, but the 
prescription I introduced in Ref. [[[(J does not. 
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For prescription II and III, we calculated A in the large-iV limit to all orders in the 
LDE, and studied its convergence properties. We found that it seems to converge to the 
correct analytic result for // > /j c , where /j c ~ 0.7. The convergence for fixed \i is extremely 
slow: the errors seem to decrease roughly like n -1 / 3 . The convergence can be accelerated 
by using the PMS criterion to choose a value of /i that depends on the order in 5. The 
PMS criterion is a polynomial equation with multiple roots, so there are various families 
of predictions for A. The family consisting of the real solutions of the PMS criterion gives 
the slowest possible convergence: the errors in A seem to decrease roughly like rT x l 2 . The 
family consisting of the complex solutions of the PMS criterion with the maximal value of 
|ImA| have nearly optimal convergence: the errors in Re A seem to decrease exponentially 
with n. Unfortunately, accurate predictions still require calculations to high orders in the 8 
expansion. To achieve 10% accuracy requires computing to about 18 th order in 5. 

For the case N = 2 relevant to Bose-Einstein condensation, we calculated A to 4 th order 
in the LDE using both prescriptions II and III. The quantity A was previously calculated 



to 4 order using prescription I in Ref. pL7 1 . The PMS criterion has complex solutions 



beginning at 3 rd order for prescription I and at 4 th order for prescriptions II and III. We 
take the values of ReA*™) for these solutions to be the predictions for A. The 3 rd and 4 th 
order predictions from prescription II are smaller than the lattice Monte Carlo result by 
about 66% and 62%, respectively. The 2 nd , 3 rd and 4 th order predictions from prescription I 
are larger than the lattice Monte Carlo result by about 131%, 85%, and 20%, respectively. 
We argue that the 20% error in the 4 th order prediction from prescription I is fortuitously 
small. Prescription I differs from prescription III by the additive term given in ([34]). The 
expansion of the factor a/1 — 8 truncated at n th order in 6 and evaluated at 5 = 1 converges 
to as n —>■ oo. The prediction from prescription I should therefore approach the prediction 
from prescription III as n increases, which suggests that it will converge to the correct answer 
from below. Because the n = 2 prediction from prescription I is greater than the correct 
result by about a factor of 2.3, there must be a crossover in n where the prediction goes from 
values greater than the correct result to smaller values. The high accuracy of the 4 th order 
result suggests that this crossover occurs at n = 5. Thus we expect higher order predictions 
from prescription I to first overshoot the correct value and then approach it from below. 

Before the lattice Monte Carlo calculation of Ref. [|, 0, there were many previous 
attempts to calculate the coefficient c in the expression ([|) for the leading order shift in T c 
for Bose-Einstein condensation. For most of these methods || ||, |10|, [11], [12|], there is no 



known systematic way of improving the calculation to make it more accurate. Such a method 
is best regarded as just a model of the critical O(N) field theory whose accuracy can only 
be determined by comparison with a more reliable calculational method. A systematically 
improvable method has two great advantages over such a model. First, the accuracy of the 
prediction can be improved with additional effort. Second, one can make a reliable estimate 
of the error in the prediction. One example of a systematically improvable method is the 
1/N expansion []3], [|, which has been used to calculate AT C to leading order and to 
next-to- leading order ||. Comparison with the lattice Monte Carlo results show that these 
first two approximations have errors of 77% and 30%, respectively. Our demonstration that 
the LDE converges in the large N limit suggests that the LDE may also be a systematically 
improvable method for N = 2. 

Although the convergence rate of the optimized LDE appears to be exponential in the 
large- N limit, it is still rather slow. One must calculate to about 18 th order in 8 to achieve 
10% accuracy. For the case N = 2 relevant to Bose-Einstein condensation, we have calculated 
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A to 4 th order in the LDE. It may be feasible to extend the calculation to 5 th order, but it 
would be very difficult to extend it to much higher order. The slow convergence in the large- 
N limit suggests that even if the LDE also converges for N = 2, a strict expansion in 5 is 
not useful for quantitative calculations. It may however be possible to use order-dependent 
mappings to change the expansion in 5 into a more rapidly converging expansion f23| . 

In conclusion, we have shown that the LDE for the quantity A that determines the 
leading order shift in T c converges in the large N limit if we use an appropriate prescription. 
If the PMS criterion is used to optimize the convergence, the errors seem to decrease expo- 
nentially in the order of the LDE. This suggests that the LDE may also be a systematically 
improvable approximation scheme in the case N = 2 relevant to Bose-Einstein condensation. 
With the use of order-dependent mappings, it may be possible to to develop the LDE into a 
general and powerful tool for quantitative calculations in superrenormalizable field theories, 
even at a critical point. 
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APPENDIX A: DIAGRAMS 



In this appendix, we give the expressions for the diagrams that contribute to A through 
order 5 4 . The diagrams for A are obtained by inserting propagator corrections into the 
diagram in Fig. [I], with the propagator corrections subtracted at momentum. In addition, 
any 2- loop self-energy subdiagram like that in Fig. ^must also be subtracted at momentum. 

We first introduce a compact notation for the massive propagator and for the functions 
of the momentum that correspond to subdiagrams in which two points are connected by 2 
or 3 lines: 

Di(p) = ~2~T — 2' ( AA ) 

D 2 (p) = [ -^— 2 - \ 2 o, (A.2) 

J q q 2 + m 2 ( p + q)^ + m 2 

D?,{p) = I 2 2 2 2 - ; ■ ~2~ ^5 (A- 3) 

Jqr q* + mr r z + rrr (p + q + rj + m z 
where / denotes the 3-dimensional integral over the momentum q : 

d3q (A.4) 



(2tt) 3 ' 

There is a simple analytic expression for the function D 2 (p): 



1 p 

D-zip) — arctan — . (A. 5) 

Anp 2m 
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The function D 3 (p) is ultraviolet divergent, but the divergence can be eliminated by a 
subtraction at momentum. The subtracted function can be evaluated analytically: 



D 3 (p)-D 3 (0) 



(47T) 



3m p 1 p + 9m 

1 arctan log — 

p 3m 2 9m 2 



(A.6) 



We proceed to give the expression for the contributions to A using the prescription II 
given in (|32"D. The contribution from the 2-loop propagator correction in Fig. |^ is 



A 2 
h 



1 



{p 



- [I>S (p) " ^3 



The contribution from the sum of the 3-loop propagator corrections in Fig. [7| is 



27 



13 = f-n&f. lDi(p + q)-D 1 (q)]D 2 2 (q) 



(p 

The contributions from the 4-loop propagator corrections in Fig. |8] are 



U 4 4 ^(^ + 2)(iV 2 + 6iV + 20) r 

4a = 8 81 

1 4 iV(iV + 2) (5iV + 22) r 
A 46 = —0 11 — I Ab 

1 
4 
1 

12 
1 



■la • 



A 4c 
A 4d 



5 i 4 N(N + 2) (5iV + 22) 



4c- 



9 

N + 2) 2 
9 



4d • 



A 4e = — <$V 
36 



Ue, 



where the momentum integrals are 



ha 

hb 

he 

hd 

he 



/ * />! (p + (g)]D 2 3 (g), 

Jp (p z ) Jq 

Ir&l D 1 (q)D 1 (q-r)D 2 (r) 
Jp {p z ) -V 

x [D x (r - p) D 2 (p-q)- D 1 (r) D 2 (q)} 
[D 1 (p + q)-D 1 (q)]D 1 (r)D 1 (q + r) 

Jp (p z ) Jqrs 

xD l (s) D l (q + s) D 2 (r - s) , 
J^J[D 2 {p + q)- D 2 {q)\ D\ (q) [D 3 (q) - D 3 



(A.7) 
(A.8) 

(A.9) 
(A.10) 

(A.ll) 
(A.12) 
(A.13) 
(A. 14) 
(A.15) 

(A.16) 

(A.17) 

(A.18) 
(A.19) 

(A.20) 
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If one uses the prescription III for A given in ([32]), the corresponding expressions for 
the diagrams are obtained by replacing each factor of 1/p 2 in the integrand by D\(p). We 
denote the integral corresponding to I n by J n . If one uses the prescription I of Ref. |j 
there is also a diagram that is th order in 5: 



A 
Jo 



NJn 



Di{p) ~ (p 



(A.21) 
(A.22) 



The 4 th in the class of diagrams that survive in the large- N limit is shown in Fig. |3|. 
The expression for the first three diagrams in the series are given in ( |A.7| ) , ( |A.9| ) and ( |A.11|) . 
The general expression for the (n + l)-loop diagram is 



2C n (N) - 



S u 



In = J, A? J, Pi (Q)}^ 1 [Di (p + q)- D l (q)} . 



(A.23) 
(A.24) 



If we use prescription III, the corresponding integral is 

J n = [ Dl(p) [ [D 2 {q)r 1 [Di{p + q)-D l {q)]. (A.25) 

J p J q 

A closed- form expression for the O(N) factor associated with this class of diagrams is 

CnW= (^)% (§)" (" -*>(" + '> , (A.26) 

where we have chosen the overall normalization such that C n (1) = 1. For n = 2,3 and 4, 
this expression reproduces the O(N) factor in (|A.7| ), ( |A.9| ) and (|A.11| ). Note that for n > 3, 
the large- limit of the coefficient is 



C n {N) 



N + 2\ n 



[1 + 0{l/N)] 



n > 3. 



(A.27) 



The corrections are suppressed by powers of 1/N. If we use the simpler expression (A^/3) rt , 
the error is of order 1/N for n N and of order A^° for n of order N. 



APPENDIX B: INTEGRALS 



In this appendix, we evaluate the momentum integrals that appear in the expressions 
for the diagrams in Appendix A. 

We first consider the integrals I n that arise if we use prescription II for A given in ([32]) . 
The 3-loop and 4-loop integrals (|A.8| ) and (|A.10|) can be evaluated analytically: 



h 
h 



6(47rfm 
1 

'24 (Att) 4 



vr 2 + 161og- + 12Li 2 (-l/3) 



(B.l) 
(B.2) 
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The remaining integrals must be evaluated numerically. Three of the 5-loop integrals can be 
reduced to integrals over a single momentum. The integral h e in ([A. 20 ) is already in such 
a form. The integrals ha in ( |A.16| ) and hd in ( |A.19| ) can be reduced to this form by first 
integrating over p : 

1 Tfl I 

w ^ + q) - Dllq)] = -- w -^, (B.3 

■^ f yD^,)-DM\ = ~^ rr (B.4 
The resulting expression for the integrals are 

/*, = ~j q D\ (?) D>(q), (B.5 

i « = -^L^ d ^- d ^w^- (B - 6 

Reducing these integrals to 1-dimensional integrals and evaluating then numerically, we 
obtain 

ha = —^6 (-0.338124209)^, (B.7 
(47r) m 6 

hd = 7-^-6 (0.0196337530) \, (B.8 
(47r) m 6 

he = -^(0.0148131532)^. (B.9 

The remaining two 5-loop integrals can be expressed as 3-loop integrals with the Mercedes- 
Benz topology. The integral hb in ( A.17|) already has this form. The integral h c in ( A. 18 
can be put in this form by integrating over p: 



m 

'4c ~ 



— / D\ (q) Di (r) D x (s) D x {s - q) D x {q - r) D 2 (r - s) . (B.10 

47T J qrs 



These integrals can be reduced to 7-dimensional integrals that are then evaluated numerically 
using the Monte Carlo method to obtain 

h b = i (-0.25) (B.ll 
(47r) m 6 

he = -^(-0.18)—. (B.12 

(4tt) 6V ; m 3 v 

The integrals J n required if A is calculated using the prescription III differ from the 
integrals J„ calculated above only in the substitution 1/p 2 — > D\{p). The 3-loop integral 
can be evaluated analytically: 

J 2 = — ^log-— . (B.13) 
2 (4tt) 3 & 4 m K J 

The 4-loop and 5-loop integrals must be evaluated numerically. Four of them can be reduced 
to integrals over a single momentum. The integrals J3 and already have this form. The 
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integral Ji a and J±d can be put into this form by first integrating over p: 

^;WP,(p + ,)- Bi (,)| = -A__£L__, ( B.i4) 

1 1 



/ D\ (p) [D 2 (p + q)-D 2 (q)} 

Jv 



32n 2 mq 



q q 

arctan arctan 

3m 2m 



(B.15) 



Reducing the final momentum integrals to a 1-dimensional integral and then evaluating it 
numerically, we obtain 

J 3 = 7A5 (-0.642143990) \, (B.16) 
(47r) m 

■ha = -^(-0.248926826)^, (B.17) 
(47r) m 6 

J M = -^(0.0165690298)^, (B.18) 



(4tt) d 

^(0.00915691397)-! 



J 4e = ^(0.00915691397)—. (B.19) 



The remaining two 5-loop integrals can be expressed as 3-loop integrals with the Mercedes- 
Benz topology. The integral in ( A. 17 ) already has this form. The integral in ([A. 18 



can be put in this form by integrating over p. These integrals can be reduced to 7-dimensional 
integrals that are then evaluated using the Monte Carlo method to obtain 

^ = ^(-0.19)^, (B.20) 
^ = ^(-0.13)^. (B.21) 



If we use the prescription I of Ref. |T7|, we must also evaluate the integral J defined in 
flOp : 

J = — -m. (B.22) 
An 

The integrals I n that survive in the large- N limit can all be reduced to integrals over a 
single momentum. After interchanging orders of integration in ( |A.24| ) and integrating over 
p, the remaining integral reduces to 

i \ n ~* r j y 2 r,/.,,M 



, dy - . [A(y)] n ~ x , (B.23) 

3tt 3 V87rmy Jo {y 2 + if 

where the function A(y) is defined in d39|). Applying the same strategy to ( |A.25|) , we obtain 
^"16^ (id L " V + 1W + 4) ^)1 (B.24, 
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TABLE I: The values of A/(— Nu/96tt 2 ) in the large- N limit at n th order in the 5 expansion using 
the prescriptions II and III for various choices of fi. Columns 1 and 4 are for fixed values of /j,. 
Columns 2 and 5 are the values at the real solutions to the PMS criterion. Columns 3 and 6 are 
the values at the complex solutions with the maximal value of |ImA|. 





II 


III 






PMS 


PMS 




PMS 


PMS 


n 


H = 1.039 


H real 


max | Im A | 


H = 0.947 


H real 


max | Im A | 


3 


0.481 


0.481 




0.456 


0.456 




4 


0.625 




0.557 ±0.156« 


0.601 




0.525 ± 0.157i 


5 


0.603 


0.611 


0.615 ±0.245i 


0.574 


0.579 


0.574 ± 0.250i 


7 


0.656 


0.675 


0.705 ± 0.339i 


0.627 


0.641 


0.653 ± 0.355i 


11 


0.709 


0.742 


0.818 ±0.400i 


0.683 


0.708 


0.755 ± 0.438i 


19 


0.761 


0.801 


0.913 ±0.416i 


0.737 


0.769 


0.851 ± 0.480i 


35 


0.809 


0.850 


0.971 ±0.410i 


0.787 


0.821 


0.917 ±0.493i 
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TABLE II: The values of A/(— u/48ir 2 ) for N = 2 using the prescriptions I, II, and III at n th order 
in the LDE with /j, evaluated at each of the solutions to the PMS criterion. 



n 


I 


II 


III 


2 


1.314 








-1.314 






3 


1.051 ±0.710i 


0.192 


0.182 




-1.920 






4 


0.679 ± 0.969i 


0.214 ± 0.084i 


0.200 ± 0.083i 




1.312 








-2.270 







26 



